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PROJECT  ABSTRACT 


By  stressing  the  sea  surface,  the  marine  atmospheric  boundary  layer 
(MABL)  can  alter  the  sea  surface  wave  field  and  so  can  produce  discernible 
signatures  on  SAR  images  of  the  ocean  (e.g.  Visecky  and  Stewart  1982).  Among 
the  resulting  signatures,  the  quasi-linear  and  cellular  microscale  patterns  still 
require  adequate  explanation.  The  ubiquitous  MABL  two-  and  three-dimensional 
convective  circulations  provide  promising  candidates  for  the  forcing  phenomena 
producing  these  signatures.  These  microscale  circulations  have  horizontal 
wavelengths  on  the  order  of  one  to  ten  times  the  boundary  layer  depth,  or 
approximately  one  to  ten  km,  and  temporal  scales  on  the  order  of  one  to  ten 
hours.  Thus,  they  produce  stress  variations  on  the  spatial  and  temporal  scales  of 
the  quasi-linear  and  cellular  SAR  signatures. 

Long  Term  Goals: 

Our  ultimate  goal  is  to  develop  methods  for  diagnosing  both  the  form  and 
the  effect  on  sea-surface  stress  patterns  of  buoyantly-forced  MABL  flows  given 
only  values  of  the  large-scale  meteorological  and  oceanographic  parameters.  As 
briefly  summarized  below,  we  continue  to  make  strong  progress  on  this  problem 
using  several  interacting,  complementary  techniques  that  range  from  data  analysis 
to  model  development. 

Near  Term  Objectives: 


Our  most  immediate  goal  is  to  determine  how  effective  the  two-and  three- 
dimensional  buoyantly  forced  MABL  circulations  with  scales  on  the  order  of  0. 1 
to  10  km  are  in  producing  microscale  patterns  of  sea-surface  stress  variability  that 
can  be  directly  linked  to  the  quasi-linear  and  cellular  SAR  signatures  that  were 
observed  during  the  Hi-Resolution  ARI  experiment  and  by  earlier,  satellite 
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mounted,  sensors.  Two  primary  objectives  must  be  met  for  us  to  achieve  these 
goals.  First,  the  surface  layer  structure  within  the  MABL  must  be  related  to  the 
sea-surface  stress  patterns  for  the  solenoidal  MABL  circulation  associated  with 
quasi-linear  sea-surface  temperature  (SST)  gradients  and  for  each  of  the  three 
microscale  atmospheric  boundary  layer  convective  forms:  the  two-dimensional 
mixed-layer  rolls,  the  three-dimensional  mixed-layer  thermals,  and  the  surface 
layer  plumes.  Second,  the  environmental  conditions  necessary  and  sufficient  for 
the  formation  of  each  of  these  buoyantly-driven  circulations  must  be  identified. 

Approach: 

The  convective  component  of  the  first  objective  has  been  completed  as  the 
MS  thesis  of  Todd  Sikora  (Sikora  1992;  Sikora  and  Young,  1993)  by  using  both 
conditional  sampling  and  composite  analysis  of  atmospheric  surface  layer 
observations  to  describe  the  occurrence,  structure,  and  sea-surface  stress  and  wind 
patterns  of  MABL  convection.  Aircraft  data  from  the  marine  stratocumulus- 
topped  boundary  layers  observed  during  Project  FIRE  [First  ISSCP  (International 
Satellite  Cloud  Climatology  Program)  Regional  Experiment]  were  used  for  the 
data  analysis.  Conditional  sampling  of  these  boundary-layer  data  identified 
updrafts  and  downdrafts,  while  composite  analysis  was  used  to  summarize  the 
spatial  (planview)  variations  within  updrafts  and  downdrafts  to  quantify  the 
average  structure  of  each.  Current  work  focuses  on  preparation  of  a  technical 
memorandum  detailing  how  to  use  these  results  as  forcing  in  ocean  wave  models. 

This  task  will  be  followed  by  a  comparison  of  these  observational  results  with 
output  from  the  new,  intermediate-order  models  described  below. 

With  the  observational  analysis  of  microscale  convection  nearing 
completion,  we  are  proceeding  to  the  next  MABL  phenomenon  of  interest  to  the 
Hi-Resolution  ARI  community,  the  MABL  solenoidal  circulation  associated  with 
the  sharp  SST  gradient  along  the  northwest  wall  of  the  Gulf  Stream. 

Observations  from  the  two  ARI  cruises  will  be  used  to  create  similar  composite 
analyses  of  this  somewhat  larger  scale,  yet  still  buoyantly  driven,  flow.  The  ship- 
based  sensors  deployed  during  both  cruises  by  Penn  State  and  collaborators  from 
Wood's  Hole  Oceanographic  Institute  will  permit  measurements  of  this  MABL 
response  to  the  ocean  front  as  well  as  the  resulting  feedback  to  the  wave  field  via 
cross-front  variations  in  the  wind  and  stress.  Cross-front  sampling  legs  and 
along-front  legs  with  two  ships  paired  across  the  front  will  yield  maximum 
information  on  the  role  of  the  SST  gradient  in  forcing  the  observed  horizontal 
variations  of  the  surface  wave  field. 

i 

Elucidation  of  the  stress  variability  expected  in  different  forcing  regimes  naumotuntat 
is  being  achieved  by  the  development  and  study  of  new  two-dimensional  and 
three-dimensional  intermediate-order  spectral  models  of  MABL  convection.  The - 
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first,  a  lower  resolution  version  following  the  boundary  layer  roll  model  of  Haack 
and  Shirer  (1992),  is  the  MS  thesis  project  of  Peter  Bromfield.  The  second  is  a 
higher  resolution  three-dimensional  convective  model  whose  development  was 
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begun  by  Julie  Schramm  and  will  be  continued  by  MS  student  Louis  Zuccarello. 
In  both  these  models,  the  boundary  conditions  are  based  on  boundary  layer 
similarity  theory  and  so  are  sufficiently  general  that  the  surface  stress  can  be 
nonzero. 

To  help  guide  comparison  of  observed  and  modeled  results,  algorithms  for 
adequately  estimating  such  chaos  measures  as  the  correlation  dimension  are 
needed  as  well  as  a  means  for  estimating  the  probable  errors  in  these  estimates. 
We  have  developed  such  an  algorithm  (Wells  et  al.  1993),  and  MS  student 
Christian  Fosmire  has  applied  this  algorithm  to  boundary  layer  wind  data. 

Recent  Tasks  Completed: 

The  conditional  sampling  and  composite  analysis  of  the  planview  patterns 
of  surface  stress  variability  caused  by  MABL  convective  updrafts  and  downdrafts 
has  been  completed  with  the  results  accepted  for  publication  (Sikora  1992;  Sikora 
and  Young  1993)  and  presented  at  the  Tenth  Symposium  on  Turbulence  and 
Diffusion  (Sikora  and  Young  1992). 

Meteorological  observations  and  infrared  satellite  images  have  been 
acquired  for  interpretation  of  an  ERS-1  SAR  feature  in  collaboration  with  Hi- 
Resolution  ARI  researcher  Robert  Beal  of  Johns  Hopkins  University's  Applied 
Physics  Laboratory  (APL). 

In  collaboration  with  Don  Thompson,  another  ARI  researcher  at  APL, 
data  acquisition  and  theoretical  explanation  have  been  completed  for  the  MABL 
front  and  solenoidal  circulation  observed  during  a  transect  of  the  Gulf  Stream 
ocean  front  by  the  RV  Vernadsky.  Mesoscale  atmospheric  simulations  of  the 
similar  flows  observed  during  the  first  Hi-Resolution  ARI  cruise  in  Fall  1991 
were  completed  in  real  time  using  the  Penn  State  mesoscale  model  MM4,  which 
was  also  used  to  provide  daily  forecast  support  to  the  cruise  planners. 

A  version  of  the  three-dimensional  intermediate-order  model  was 
integrated  and  preliminary  results  showed  that  the  forms  of  the  stress  and  heat 
flux  profiles  were  correct  and  matched  those  observed  in  the  stratocumulus- 
topped  MABL  by  Briimmer  and  Busack  (1990).  The  energy  budget  properties 
were  not  completely  correct,  however,  and  appear  to  be  caused  in  part  by  an 
improper  handling  of  the  pressure  gradient  term  in  the  equation  of  motion. 
Because  this  term  is  much  easier  to  handle  in  a  two-dimensional  model,  we 
undertook  the  development  of  such  a  model  as  the  MS  thesis  project  for  Peter 
Bromfield.  We  also  began  a  reexamination  of  the  formulation  of  the  three- 
dimensional  model;  the  corrected  formulation  will  be  implemented  by  MS  student 
Louis  Zuccarello. 

The  correiution  dimension  algorithm  of  Wells  et  al.  (1993)  was 
supplemented  with  means  for  estimating  the  relative  error  in  the  estimate  of  this 
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dimension.  Christian  Fosmire  applied  this  algorithm  to  boundary  layer  wind  data 
and  determined  that,  when  an  estimate  of  the  relative  error  is  included,  the 
dimension  values  given  by  either  horizontal  wind  component  are  the  same  as 
those  expected  from  theory;  without  considering  these  error  estimates,  the 
dimension  values  would  appear  to  be  different  (Fosmire  1993). 

Results: 

The  published  analyses  of  the  composite  structure  of  MABL  updrafts  and 
downdrafts  (Sikora  1992;  Sikora  and  Young  1992,  1993)  provide  the  first-ever 
quantitative  description  of  the  planview  patterns  of  air/sea  fluxes  of  momentum, 
heat,  and  moisture.  Figure  1  shows  the  corresponding  composite  planview 
perturbation  wind  field  for  MABL  convective  downdrafts,  the  phenomenon  that 
drives  "cat's  paw"  wave  patterns  on  the  sea  surface.  These  results  can  be  used  by 
ocean-wave  modelers  to  quantify  the  horizontal  variability  of  atmospheric  forcing 
on  spatial  scales  on  the  order  of  0.1  to  1.0  km. 

Preliminary  theoretical  analysis  suggests  that  a  solenoidal  circulation,  akin 
to  the  classic  sea-breeze  circulation,  causes  the  significant  MABL  wind  variations 
observed  in  the  vicinity  of  the  northwest  wall  of  the  Gulf  Stream.  This  common 
phenomenon  is  revealed  by  a  cloud  line  in  meteorological  satellite  and  Hi- 
Resolution  cruise  observations  as  well  as  by  the  band  of  altered  wind  direction  in 
the  MM4  mesoscale  model  forecasts  and  the  RV  Vernadsky  transect.  Further 
analysis  of  observations  from  the  second  Hi-Resolution  ARI  cruise  and 
quantitative  diagnostics  of  the  MM4  results  will  be  required  to  complete  the 
description  of  this  phenomenon,  its  effect  on  the  sea-surface  stress  field,  and  the 
range  of  conditions  under  which  it  can  occur. 

A  manuscript.  Wells  et  al.  (1993),  summarizing  the  new  algorithm  for 
estimating  the  correlation  dimension  and  its  probable  error  is  nearly  completed 
and  will  be  submitted  to  Physica  D  in  1993.  A  preliminary  description  of  the 
results  in  this  paper  was  presented  at  the  Spring  1992  AGU  conference  in 
Montreal,  Canada  (Wells  et  al  1992).  This  algorithm  has  been  applied  to 
boundary  layer  wind  observations  by  Fosmire  (1993)  and  is  ready  to  be  applied  to 
other  modeled  and  observed  MABL  datasets. 

Accomplishments: 

A  graduate  student  (Todd  Sikora)  has  completed  his  MS  thesis  (Sikora, 
1992)  on  the  patterns  of  horizontal  variability  in  sea-surface  forcing  caused  by 
MABL  convection.  These  results  have  been  presented  by  George  Young  at  the 
April  14-16,  1992  Hi-Resolution  Workshop  at  APL,  by  Todd  Sikora  at  the  Tenth 
Symposium  on  Turbulence  and  Diffusion  in  Portland,  Oregon,  September  29- 
October  2,  1992  (Sikora  and  Young  1992),  and  have  been  accepted  for 
publication  in  Boundary-Layer  Meteorology  (Sikora  and  Young  1993).  Todd 
Sikora  and  George  Young  have  also  pursued  a  meteorological  explanation  for  an 
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ERS-1  SAR  feature  noted  by  Robert  Beal.  A  new  graduate  student  (Sean 
Sublette)  is  being  taught  to  analyze  mesoscale  model  output  in  conjunction  with 
Hi-Resolution  cruise  observations  to  develop  a  composite  description  of  the  flow 
and  stress  fields  associated  with  the  MABL  solenoidal  circulation  that  is  forced 
by  the  Gulf  Stream  ocean  front.  Theoretical  insight  into  this  frequently  occurring 
phenomenon  has  been  provided  to  interested  ARI  investigators,  Don  Thompson 
and  Farid  Askari.  Also,  George  Young  helped  plan  the  next  Hi-Resolution  cruise 
at  the  October  27-29,  1992  workshop  at  NI^.  Graduate  students  Peter  Bromfield 
and  Louis  Zuccarello  are  developing  two-  and  three-dimensional  models  of  the 
stress  variability  caused  by  MABL  circulations;  these  models  are  extensions  of 
the  first  version  coded  by  research  associate  Julie  Schramm.  Preliminary  results 
from  this  version  were  presented  by  Nels  Shirer  at  the  April  14-16,  1992  Hi- 
Resolution  Workshop.  Finally,  MS  student  Christian  Fosmire  has  applied  the 
Wells  et  al.  (1993)  algorithm  to  boundary  layer  wind  data;  he  is  writing  his  thesis 
(Fosmire  1993)  for  an  expected  May  1993  graduation.  This  algorithm  was  also 
presented  by  Christian  Fosmire  at  ^e  Spring  1992  AGU  meeting  in  Montreal, 
Canada  (Wells  era/.  1992). 
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Figure  1  Planview  of  composite  horizontal  wind  in  MABL  convective 
downdrafts  in  the  surface  layer. 
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Abstract.  Air/sea  flux  variability  on  horizontal  scales  from 
50  m  to  several  km  results,  in  part,  from  the  presence  of 
coherent  convective  structures  within  the  atmospheric  boundary 
layer.  The  horizontal  distribution  of  fluxes  within  these 
convective  updrafts  and  downdrafts  is,  therefore,  central  to 
studies  of  air/ sea  interaction  and  remote  sensing  of  sea 
surface  wind  and  wave  fields.  This  study  derives  these  flux 
patterns  from  observations  of  the  Marine  Atmospheric  Surface 
Layer  (MASL) . 

Research  aircraft  flights  through  the  MASL  provide  an 
optimal  means  for  sampling  large  numbers  of  the  above- 
mentioned  coherent  structures.  The  NCAR  Electra  flew  numerous 
legs  through  the  MASL  at  a  height  of  50  m  during  the  1987 
stratocumulus  phase  of  Project  FIRE  (First  ISSCP 
(International  Satellite  Cloud  Climatology  Program)  Regional 
Experiment) .  In  situ  measurements  from  these  legs  serve  as 
the  dataset  for  this  paper.  The  data  are  processed  in  such  a 
way  as  to  retain  only  the  turbulence  fluctuations. 
Conditional  sampling,  based  on  the  vertical  velocity  field, 
results  in  the  isolation  of  convective  updrafts  and 
downdrafts.  Compositing  of  the  data  for  these  two  classes  of 
convective  drafts  results  in  horizontal  planviews  of  the 
vertical  fluxes  of  buoyancy,  absolute  humidity,  along-mean- 
wind  component  of  momentum,  and  vertical  velocity.  To  ensure 
dynamical  similarity,  these  horizontal  planviews  are  oriented 
in  a  coordinate  system  aligned  with  the  mean  wind. 
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1.  Introduction 

The  marine  atmospheric  surface  layer  (MASL)  (Z  <  0.1 
where  is  the  boundary  layer  depth)  (Stull,  1989)  plays  a 
major  role  in  air/ sea  interaction.  It  is  the  MASL  that 
couples  the  sea  surface  to  the  overlying  atmospheric  boundary 
layer  through  vertical  eddy  fluxes  (hereafter  referred  to 
simply  as  fluxes) ,  such  as  those  of  buoyancy,  moisture,  and 
momentum.  While  these  fluxes  are  present  in  both  stable  and 
unstable  environments,  they  tend  to  be  enhanced  by  convective 
elements,  as  is  implied  by  bulk  aerodynamic  parameterization 
(e.g.  Liu  et  al.,  1979).  For  this  reason,  much  of  the  recent 
research  in  this  area  has  concentrated  on  the  convective 
marine  atmospheric  surface  layer  (CMASL) . 

CMASL  fluxes  are  realized  through  sub-mesoscale 
convective  updrafts  (CUs)  and  downdrafts  (CDs) .  These 
features  have  diameters  on  the  order  of  tens  to  thousands  of 
meters  (Lenschow  and  Stephens,  1980) .  Khalsa  and  Greenhut 
(1985)  found  that  within  a  central  Pacific  CMASL,  such 
features  were  responsible  for  75%  of  the  total  flux  of  heat, 
moisture,  and  momentum.  This  importance  provides  the 
motivation  to  explore  further  the  flux  characteristics  of 
CMASL  drafts. 

Various  methods  may  be  employed  to  study  CMASL  drafts. 
Qualitatively,  this  can  be  accomplished  by  simply  observing 
the  effects  their  associated  fluxes  have  on  the  environment. 
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On  a  human  scale,  anyone  who  has  enjoyed  a  day  of  sailing, 
seen  wisps  of  sea  fog,  or  witnessed  cats  paws  rippling  across 
a  body  of  water  has  personally  sensed  the  impact  of  CMASL 
drafts  on  the  environment.  Another  example  of  this 
manifestation  can  be  seen  in  Figures  l  and  2.  Figure  1  is 


PsO-V  ClOcx'A 


Fig.  1.  Visible  image  from  Fig.  2.  X-band  real 
the  Kosmos  1500  satellite.  aperture  radar  image  of  the 

underlying  sea  surface. 

a  visible  image  from  the  Kosmos  1500  (Okean)  satellite.  The 
image  is  centered  at  26°  N  125.5°  W  and  is  dated  July  ll, 
1984.  It  shows  common  boundary  layer  cumuliform  clouds, 
horizontal  scales  of  which  are  on  order  of  the  boundary  layer 
depth,  induced  by  fluxes  of  heat  and  moisture  within  CMASL- 
born  CUs.  Figure  2  is  the  corresponding  X-band  real  aperture 
radar  image  of  the  underlying  sea  surface.  It  depicts  a 
perturbed  sea  surface  wave  pattern  resulting  from  CMASL 
momentum  flux  patterns  driven  by  flow  into  CUs  and  out  of  CDs. 

In  order  to  obtain  data  on  the  effects  and  structure  of 
CUs  and  CDs  in  a  more  quantitative  sense,  various  techniques 
have  been  employed,  including  observations  from  a  sensor 
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equipped  catamaran  (Dorman  and  Mollo-Christensen,  1973)  ,  the 
use  of  acoustic  sounders  (Gaynor  and  Mandics,  1978) ,  and 
observations  of  sea  gulls  in  flight  (Woodcock,  1940,  1975) . 
Data  for  the  quantitative  study  in  this  paper  are  obtained 
from  NCAR  Electra  aircraft  flights  during  the  First  ISCCP 
(International  Satellite  Cloud  Climatology  Project)  Regional 
Experiment  (FIRE) .  Flights  were  conducted  over  a  several 
hundred  kilometer  region  west  of  the  southern  California  coast 
during  June  and  July,  1987.  For  detailed  reviews  of  the  FIRE 
project,  see  Albrecht  et  al.  (1988)  and  Kloesel  et  al.  (1988). 

The  traditional  way  to  investigate  the  flux 
characteristics  of  boundary  layer  features  is  through 
conditional  sampling  and  composite  analysis,  examples  of  which 
are  discussed  below.  Simply  put,  conditional  sampling  uses  an 
indicator  function  to  isolate  from  a  dataset  features  of 
interest  (e.g.  CUs  and  CDs).  Composite  analysis  then  combines 
information  from  these  selected  features  to  determine  their 
average  spatial  structure. 

In  the  past,  various  studies  have  gathered  valuable 
information  concerning  boundary  layer  features  using  these 
sampling  and  analysis  techniques.  Wilczak  and  Tillman  (1980) 
and  Wilczak  (1984)  used  Boulder  Atmospheric  Observatory  tower 
data  to  produce  detailed,  land-based,  statistics  for 
convective  atmospheric  surface  layer  features  known  as 
temperature  ramps.  Young  (1988)  also  investigated  land-based 
turbulence  at  the  Boulder  Atmospheric  Observatory.  He 
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conditionally  sampled  aircraft  data,  based  upon  vertical 
velocity  perturbations  and  mixed  layer  spectra  of  vertical 
velocity  and  temperature,  in  order  to  analyze  thermals  within 
the  convective  boundary  layer.  Lenschow  and  Stephens  (1980) 
used  NCAR  Electra  aircraft  data  to  generate  information  on  the 
structure  of  marine  atmospheric  boundary  layer  (MABL)  thermals 
using  humidity  as  an  indicator  function.  Greenhut  and  Khalsa 
(1982)  and  Khalsa  and  Greenhut  (1985)  also  used  aircraft  data 
to  study  properties  of  drafts  in  the  MABL  but,  as  in  Young 
(1988)  ,  they  used  vertical  velocity  as  an  indicator  of  a 
feature.  In  addition  to  observational  studies,  large-eddy 
simulation  models,  such  as  that  of  Schumann  and  Moeng  (1991)  , 
have  been  used  to  provide  boundary  layer  datasets  for 
conditional  sampling  and  compositing. 

While  previous  work  has  resulted  in  detailed  vertical 
profiles  of  properties  of  MABL  updrafts  and  downdrafts, 
information  on  planview  flux  variability  of  CMASL  CUs  and  CDs 
is  lacking.  To  fill  this  need,  this  study  utilizes 
conditionally  sampled  and  composited  eddy  correlated  FIRE 
aircraft  data  to  investigate  the  typical  CMASL  CU  and  CD 
planview  patterns  of  vertical  velocity  flux,  w'w',  buoyancy 
flux,  w'T^',  absolute  humidity  flux,  w'r',  and  along-mean-wind 
momentum  flux,  w'u' .  Note  that  while  CUs  and  CDs  are 
segregated  in  this  study,  the  techniques  used  for  their 
analyses  are  identical. 
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2 .  Procedures 

The  use  of  NCAR  Electra  turbulence  data  from  project  FIRE 
is  advantageous  for  the  type  of  research  discussed  here 
because  a  properly  instrumented  aircraft  provides  a  means  for 
collecting  high  resolution  turbulence  data  from  a  large  number 
of  CMASL  drafts.  While  each  draft  can  be  penetrated  only 
once,  penetrations  of  many  different  drafts  can  be  made  from 
numerous  angles.  In  addition,  the  turbulence  data  from  a 
draft  sampled  by  an  aircraft  can  be  looked  upon  as  a  snapshot 
of  that  feature  (see  sub-section  2.3.).  The  result,  after 
data  processing,  conditional  sampling,  and  compositing,  is  a 
statistically  robust  planview  flux  description  of  a  typical 
CMASL  draft. 

2.1.  FLIGHT  INFORMATION 

Of  the  numerous  available  datasets  collected  from  Electra 
flight  legs  during  FIRE,  only  the  20  used  in  this  study  were 
conducted  within  the  CMASL.  These  legs  averaged  50  km  (500 
sec)  in  length  (duration) .  These  legs  were  flown  on  seven 
different  days,  during  which  the  average  boundary  layer  depth 
was  920  m.  The  legs  were  flown  at  a  height  of  50  m,  well 
within  the  depth  of  the  surface  layer.  The  ratio  of  the 
boundary  layer  depths  to  the  Monin-Obukhov  lengths  ranged  from 
-0.62  to  -37.56,  indicative  of  slightly  to  moderately  unstable 


6 


MABLs.  While  flux  intensity  differences  for  CUs  and  CDs  are 
seen  within  this  stability  range,  the  "flux-shapes”  (the 
patterns  of  spatial  variation  of  flux  magnitude  across  a 
draft)  are  similar.  This  pattern  similarity  in  the  presence 
of  intensity  differences  is  exactly  analogous  the  geometric 
similarity  of  "similar  triangles"  of  varying  sizes. 
Structural  differences  are,  however,  observed  for  stable  and 
extremely  unstable  cases.  Insufficient  data  prevented  the 
extension  of  the  composite  analysis  to  these  other  stability 
ranges . 

2.2.  INSTRUMENTATION  AND  DATA  PROCESSING 

The  Electra  CMASL  20  hz  turbulence  data  used  in  the 
present  study  include  vertical  velocity,  w,  air  temperature, 
T,  absolute  humidity,  r,  and  the  aircraft-oriented  components 
of  the  wind,  u,  and,  v.  This  sampling  rate  combined  with  the 
100  m/s  aircraft  speed  yields  a  minimum  resolvable  wavelength 
of  10  m.  A  review  of  the  aircraft  instrumentation  of 
significance  to  this  study  can  be  found  in  Nucciarone  and 
Young  (1991).  Data  processing  for  mean,  trend,  and  spike 
removal,  as  well  the  derivation  of  buoyancy  flux,  parallels 
that  of  Moyer  and  Young  (1991).  In  order  to  eliminate 
contributions  by  mesoscale  phenomena  in  the  data,  high-pass 
filtering  is  employed  following  the  techniques  of  Young 
(1987,1988) . 
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An  example  of  the  spectral  response  of  the  filter  can  be 
seen  in  Figure  3.  The  figure  shows,  for  a  cosine  wave,  how 


Fig.  3.  An  example  of  the  spectral  response  of  the  filter. 

the  ratio  of  the  wave's  filtered  variance  to  its  unfiltered 
variance  varies  with  the  ratio  of  wavelength  to  Z^.  In  this 
example,  the  cutoff  wavelength  chosen  is  0.1  Z^.  At  this 
wavelength,  50%  of  the  wave's  original  variance  remains.  At 
longer  wavelengths,  the  amount  of  the  wave's  original  variance 
remaining  after  filtering  diminishes  rapidly. 

FIRE  spectra  of  Nucciarone  and  Young  (1991)  indicate  that 
1.5  Z^  is  the  appropriate  mesoscale  subrange  cutoff  wavelength 
for  use  in  the  filter.  The  net  result  of  this  data 
processing,  for  any  one  variable,  is  a  sub-mesoscale 
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perturbation  series  (denoted  by  primed  variables)  suitable  for 
computing  eddy  correlation  statistics. 

Coordinate  rotations  are  performed  on  the  perturbation 
wind  components  yielding  one  component  in  the  direction  of  the 
mean  wind,  u' ,  and  another  aligned  with  positive  values  90°  to 
the  left  of  the  mean  wind,  v' .  Unlike  that  seen  in  Geernaert 
and  Hansan  (1992) ,  where  orographic  distortion  of  the  low 
level  flow  may  be  occurring,  the  cross-mean-wind  components  of 
the  momentum  fluxes  of  this  study  are  both  statistically  and 
physically  insignificant  when  compared  with  the  a long -mean- 
wind  component  (i.e.  mean  <  standard  deviation  for  w'v'  and 
w'v'  «  w'u').  For  this  reason,  the  cross-mean-wind  momentum 
flux  will  not  be  discussed. 

2.3.  CONDITIONAL  SAMPLING 

In  order  to  distinguish  updrafts  and  downdrafts,  event 
criteria  based  on  perturbation  vertical  velocity  are  required. 
Greenhut  and  Khalsa  (1982)  and  Young  (1988)  show  that  w'  event 
criteria,  incorporating  germane  horizontal  scales  and 
magnitude  thresholds,  allow  for  the  proper  detection  of 
boundary  layer  features  of  interest. 

In  order  to  distinguish  CUs  and  CDs,  the  w'  series  is 
first  band-pass  filtered  to  retain  only  the  convective  scales 
(Young  1988)  .  The  band-pass  filter  is  a  variation  on  that 
discussed  in  sub-section  2.2.,  combining  both  a  high-pass  and 
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low-pass  stage.  It  is  designed  to  eliminate  both  the 
mesoscale  and  inertial  subrange  contributions  to  the  w' 
series,  while  still  preserving  the  CUs  and  CDs  of  the  energy 
containing  subrange.  As  discussed  in  sub-section  2.2.,  1.5 
is  used  as  the  mesoscale  subrange  cutoff  wavelength.  The 
choice  of  an  inertial  subrange  cutoff  wavelength  is  0.1  Z^. 
As  in  Young  (1988) ,  this  choice  eliminates  scales  of  motion 
associated  with  the  inertial  cascade. 

Given  the  average  Zj^  of  this  study,  a  minimum  draft  width 
of  50  m  is  chosen  so  that  the  narrowest  acceptable  event 
(CU/CD  couplet)  corresponds  with  the  0.1  Z^  inertial  subrange 
cutoff  wavelength  used  in  the  band-pass  filter.  All  drafts 
not  meeting  the  minimum  width  criteria  are  rejected  because  of 
their  inertial  subrange  character. 

The  literature  mentioned  in  section  1.  provides  numerous 
examples  of  conditional  sampling  magnitude  thresholds  designed 
to  distinguish  features  of  interest  from  their  environment. 
Inspection  of  the  FIRE  w'  data  series  suggests  that  a  v' 
threshold  of  ±  0.1  m/s  eliminates  disorganized  areas  of  weak 
ascent  or  decent.  Any  data  points  for  which  w'  <0.1  m/s  are 
therefore  considered  not  to  be  part  of  significant  convective 
features  and  are  so  rejected. 

CUs  (CDs)  are  then  defined  to  be  features  within  the 
band-pass  filtered  w'  data  series  meeting  these  minimum  width 
and  magnitude  criteria.  The  total  number  of  CUs  sampled  and 
used  in  compositing  is  2839,  averaging  98  m  in  width  and 
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occupying  28%  of  the  data  series.  The  total  number  of  CDs 
sampled  and  used  in  compositing  is  3062,  averaging  107  m  in 
width  and  occupying  33%  of  the  data  series.  Given  the  above 
widths  and  an  aircraft  velocity  of  100  m/s,  it  can  be  shown 
that  Taylor's  Hypothesis  is  valid  (Stull  1989) .  Draft 
composites  of  eddy  correlation  statistics  can  therefore  be 
looked  upon  as  typical  planview  flux  snapshots  of  a  CU  or  a 
CD. 


2.4.  COMPOSITING  OF  EDDY  CORRELATION  STATISTICS 

After  the  isolation  of  CUs  and  CDs,  eddy  correlation 
fluxes  of  vertical  velocity,  buoyancy,  absolute  humidity,  and 
along-mean-wind  momentum  are  calculated  from  the  high-pass 
filtered  dataset,  for  each  draft.  Note  that,  unlike  the  band- 
passed  series  used  to  locate  the  sub-mesoscale  CUs  and  CDs, 
this  series  retains  the  inertial  subrange  contributions.  Only 
the  mesoscale  contributions  are  filtered  out,  as  discussed  in 
sub-section  2.2..  In  order  to  develop  typical  planview  flux 
snapshots,  the  fluxes  from  all  the  drafts  are  composited 
together  using  the  series  of  averaging  calculations  and 
grouping  techniques  discussed  below. 

Each  draft  is  first  divided  into  3  bins  of  equal  length: 
the  aircraft  entry  region  of  the  draft,  bin  (a)  ,  the  middle  of 
the  draft,  bin  (b)  ,  and  the  aircraft  exit  region  of  the  draft, 
bin  (c)  .  The  average  CU  and  CD  flux  patterns  for  each  leg  are 
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then  found  by  calculating  bin  average  eddy  correlation 
statistics  as  follows.  The  flux  of  any  variable  y  by  variable 
w'  is  found  by  first  taking  the  sum  of  their  product  over  all 
points  of  all  like  bins  [bin  (a) ,  (b) ,  or  (c) ]  over  all  like 
drafts  (CUs  or  CDs)  within  the  leg.  This  sum  is  then  divided 
by  the  number  of  data  points  in  all  like  bins  of  all  like 
drafts  within  the  leg. 

The  resulting  40  leg  average  drafts  (LADs)  from  the  20 
flight  legs  are  then  separated  into  one  of  3  equal  angular- 
width  groups.  The  grouping  of  an  LAD  depends  on  the  path 
taken  by  the  aircraft,  through  the  draft,  relative  to  the  mean 
wind.  In  polar  coordinates  with  the  mean  wind  vector  aligned 
towards  0°,  these  relative  heading  groups  are:  along-mean- 
wind  (330°  to  30°  and  150°  to  210°),  cross-mean-wind  (60°  to 
120°  and  240°  to  300°)  ,  and  diagonal-to-the-mean-wind  (the 
remaining  relative  heading  ranges) . 

To  maintain  wind-relative  consistency  in  bin  naming,  flux 
data  for  bins  (a)  and  (c)  are  switched  for  those  drafts  of  the 
along-mean-wind  and  diagonal  groups  whose  relative  headings 
oppose  the  mean  wind.  The  same  is  done  for  LADs  of  the  cross¬ 
mean-wind  group  whose  relative  heading  range  is  aligned  to  the 
left  of  the  mean  wind.  Note  that  while  there  are  twice  as 
many  relative  heading  ranges  for  the  diagonal  group  as  for  the 
other  two,  the  angular  width  of  each  heading  range  in  the 
diagonal  group  is  only  half  as  large  as  for  the  other  two. 
Thus,  the  angular  coverage  of  all  3  groups  is  equal. 


12 


The  two  remaining  relative-heading  ranges  of  the  diagonal 
group  (30°  to  60°  and  300°  to  330°)  can  be  merged  using 
symmetry  arguments.  Wilczak  (1984)  shows  mirror  symmetry  of 
the  horizontal  perturbation  wind  field  about  an  axis  aligned 
with  the  mean  wind  extending  through  the  center  of  his  large- 
scale  eddy  (LSE  which  is  a  CU/CD  couplet) .  A  similar  symmetry 
can  be  seen  in  the  current  dataset  by  examination  of  the 
fluxes  of  the  cross-mean-wind  group  composites  found  in 
section  3..  Comparison  of  the  data  from  the  two  relative 
heading  ranges  of  the  diagonal  group  (not  shown)  demonstrates 
symmetry  across  the  axis  of  the  mean  wind  for  both  the  upwind 
and  downwind  ends  of  the  drafts.  Taking  these  symmetries  into 
account,  flux  data  of  like  LAD  dow.^xnd  bins  of  the  diagonal 
group  are  composited  together  as  are  flux  data  of  upwind  bins. 
Figure  3  shows  the  spitial  relation  of  the  three  groups 
withxn  a  planview  draft.  Bins  are  labeled  (a) ,  (b) ,  or  (c) . 
Upper  case  letters  indicate  the  group  to  which  each  bin 
belongs;  the  a long-mean-wind  group  (A) ,  the  cross-mean-wind 
group  (C) ,  and  the  diagonal-to-the-mean-wind  group  (D) .  The 
mean  wind  vector  is  assumed  to  be  directed  towards  the  top  of 
the  figure. 

The  process  of  compositing  flux  data  of  like  LAD  bins 
into  a  flux  for  the  corresponding  bins  of  a  group  average 
draft  (GAD)  will  now  be  explained.  First,  because  the  center 
bin  for  all  groups  corresponds  to  the  same  part  of  the  draft, 
the  average  bin  (b)  flux  over  all  like  LADs  over  all  groups 
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Fig.  4.  The  spatial  relation  of  the  three  groups  within  a 

planview  draft. 


(bin  (b)  average)  may  be  used  to  represent  the  flux  in  the 
center  bin  of  all  groups.  As  stated  in  sub-section  2.1., 
drafts  of  each  leg  are  structurally  similar  in  terms  of  their 
flux-shape.  However,  there  is  a  need  to  correct  for  mean 
draft  intensity  differences  between  like  bins  of  different 
LADS.  It  follows  then  that  LAD  bin  (a)  and  (c)  fluxes  are 
rescaled  before  any  averaging  to  produce  GADs.  This  rescaling 
uses  a  correction  coefficient  based  on  the  need  for  inter-leg 
similarity  in  bin  (b)  ,  described  above.  This  correction 
coefficient  is  bin  (b)  average  divided  by  the  respective  LAD 
bin  (b)  .  This  approach  is  robust  because  the  LAD  bin  (b) 
fluxes  of  all  the  quantities  studied  are  generally  large  and 
of  same  sign. 

The  final  planview  composite  is  found  by  simply 
overlaying  the  along-mean-wind,  cross-mean-wind,  and  diagonal- 


to-the-mean  wind  (again  using  symmetry  for  the 
diagonal)  GADs  for  each  statistic  discussed. 
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3 .  Observational  Results 

Planview  draft  flux  composites  of  vertical  velocity, 
buoyancy,  absolute  humidity,  and  the  along-mean-wind  component 
of  momentum,  along  with  corresponding  bin  standard  deviations, 
are  presented  in  Tables  I  through  VIII.  Each  table  is  a  three 
by  three  array  of  group  bin  fluxes.  Bin  location  in  this 
array  corresponds  to  its  position  within  the  composite  draft. 
In  all  tables,  the  mean  wind  is  directed  towards  the  top  of 
the  table.  For  each  composite,  then,  the  along-mean-wind 
group  axis  is  aligned  from  the  upwind  center  to  downwind 
center,  the  cross-mean-wind  group  axis  is  aligned  from  the 
middle  left  to  middle  right,  and  the  diagonal  group  makes  up 
the  remainder. 

Note  that  the  middle  center  bin  of  all  tables  is  its 
respective  bin  (b)  average  of  sub-section  2.4..  Recall  also 
that  the  middle  center  bin  is  used  to  rescale  all  other  bins 
within  a  flux  composite.  For  these  reasons,  the  majority  of 
the  uncertainty  associated  with  any  flux-shape  and  intensity 
related  inter leg  differences,  associated  with  inter leg 
differences  in  the  middle  center  bins,  is  forced  into  the 
perimeter  of  the  composite  array.  It  follows  then  that 
standard  deviations  for  the  middle  center  bins  are  not  given. 


3.1.  VERTICAL  VELOCITY  FLUX 


Vertical  velocity  dictates,  to  a  large  extent,  the  other 
flux  patterns  discussed  in  this  study.  Not  only  is  it  used 
for  the  defining  of  events,  but  also  it  is  the  flux  variable 
doing  the  work,  so  to  speak,  by  advecting  the  other 
quantities.  The  planview  CU  w'w'  composite  is  presented  in 
Table  I. 

TABLE  I 

CU  vertical  velocity  flux  composite  and 
corresponding  standard  deviations. 


Because  all  bin  standard  deviations  are  much  smaller  than 
the  bin  fluxes,  there  is  statistical  confidence  in  the  entire 
CU  w'w'  composite.  The  strongest  flux  region  is  found  along 
the  along-mean-wind  group  axis  with  the  largest  magnitude  of 
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flux  being  at  the  center  of  the  composite.  This  shape 
compares  well  with  that  found  at  50  m  over  land  by  Wilczak 
(1984) .  Examining  the  perimeter  of  the  composite,  the  flux 
magnitude  along  the  upwind  side  of  the  draft  tends  to  be 
larger  than  that  along  the  downwind  side,  although  these 
asymmetries  are  less  statistically  significant  than  the  radial 
variation.  The  weakest  flux  is  found  in  the  middle  right  and 
left  bins. 

The  planview  CD  vertical  velocity  flux  composite  is 
shown  in  Table  II.  As  with  the  CU  composite,  there  is 


TABLE  II 

CD  vertical  velocity  flux  composite  and 
corresponding  standard  deviations. 


w'w'  (m^s“^) 
(std.  dev.) 

Left 

Center 

Right 

Downwind 

0.160 

(0.011) 

0.161 

(0.018) 

0.160 

(0.011) 

Middle 

0.171 

(0.024) 

0.204 

0.173 

(0.024) 

Upwind 

0.162 

(0.009) 

0.166 

(0.019) 

0.162 

(0.009) 
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statistical  confidence  in  all  bin  fluxes.  While  the  largest 
magnitude  of  w'w'  is  again  found  in  the  center  middle  bin,  the 
over  all  flux  composite  is  weaker  than  that  of  the  CU.  This 
finding  is  in  agreement  with  those  of  Greenhut  and  Khalsa 
(1982)  for  their  coherent  downdrafts.  As  with  the  CU 
composite,  other  aspects  of  the  w'w'  pattern  are  less 
statistically  significant.  The  strongest  flux  region  is  found 
along  the  cross-mean-wind  group  axis.  Flux  along  the  upwind 
side  of  the  draft  is  slightly  larger  than  that  along  the 
downwind  side. 

3.2.  BUOYANCY  FLUX 

The  buoyancy  flux  of  the  current  study  has  the  same  sign 
and  approximates  to  a  large  extent  the  magnitude  of  the  heat- 
flux  w'T'.  The  planview  CU  buoyancy  flux  composite  is 
presented  in  Table  III.  While  the  composite  corners  lack 
statistical  significance,  there  is  statistical  confidence  in 
that  part  of  the  composite  consisting  of  the  along-mean-wind 
and  cross-mean-wind  groups.  Along  these  axes,  general 
symmetries  exist  with  strongest  w'T^'  being  at  the  middle 
center  of  the  composite.  All  fluxes  within  the  composite  are 
down-gradient.  This  finding  is  in  agreement  with  those  of 
Khalsa  and  Greenhut  (1985)  for  their  lowest  levels. 

The  planview  CD  buoyancy  flux  composite  is  presented  in 
Table  IV.  Statistical  significance  in  the  w'T^'  pattern  is 
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similar  to  that  of  the  CU.  The  along-mean-wind  and  cross¬ 
mean-wind  group  axes  are  more  symmetric  than  in  the  CU  while 
the  diagonal-to-the-mean-wind  group  axis  is  even  more 


TABLE  III 

CU  buoyancy  flux  composite  and  corresponding 
standard  deviations. 


w'Ty'  (ms‘^°C) 
(std.  dev.) 

Left 

Center 

Right 

Downwind 

0.024 

(0.021) 

0.010 

(0.002) 

0.024 

(0.021) 

Middle 

0.010 

(0.002) 

0.018 

0.011 

(0.002) 

Upwind 

0.006 

(0.008) 

0.014 

(0.002) 

0.006 

(0.008) 

asymmetric  in  the  along-mean-wind  direction  than  in  the  CU. 

The  pattern  of  the  CD  w'T^'  magnitude  is  similar  to  that 
of  the  CU.  However,  not  all  fluxes  are  positive.  Fluxes 
found  within  the  upwind  left  and  right  bins,  although  not 
statistically  significant,  are  counter-gradient.  In  contrast 
and  as  will  be  seen  in  sub-section  3.3.,  all  fluxes  composing 
the  CD  absolute  humidity  composite  CU  of  this  study  are  down- 
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gradient.  Keeping  in  mind  the  lack  of  statistical  confidence, 
the  upwind  left  and  right  bins  of  for  the  CD  composite, 
then,  depart  from  findings  of  Khalsa  and  Greenhut  (1985)  that 
warm  dry  downdrafts  have  negative  buoyancy  flux  only  at  levels 


TABLE  IV 

CD  buoyancy  flux  composite  and  corresponding 
standard  deviations. 


w'Ty"  (ms"^°C) 

( std .  dev . ) 

Left 

Center 

Right 

Downwind 

0.009 

(0.005) 

0.005 

(0.001) 

0.009 

(0.005) 

Middle 

0.004 

(0.001) 

0.006 

0.004 

(0.001) 

Upwind 

-0.002 

(0.009) 

0.004 

(0.001) 

-0.002 

(0.009) 

higher  than  the  surface  layer.  Similar  to  what  is  suggested 
in  Khalsa  and  Greenhut  (1985)  for  some  of  their  anomalous 
patterns,  the  counter-gradient  w'T^'  of  this  study  may  result 
from  entrainment  of  buoyant  CU  air  into  the  upwind  side  of  the 
CD.  The  positive  buoyancy  flux  portion  of  the  CD  w'T^' 
composite  of  this  study  compares  well  with  findings  in  the 
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lowest  levels  of  Khalsa  and  Greenhut  (1985). 

Of  final  note:  a  surface  layer  counter-gradient  heat 

flux  is  found  in  the  upwind  region  of  the  ensemble  downdraft 
of  Wilczak's  (1984)  vertical  cross  section.  Thus,  the 
finding  of  Wilczak  (1984)  compares  well  with  the  current 
study . 

3.3.  ABSOLUTE  HUMIDITY  FLUX 


The  planview  CU  absolute  humidity  flux  composite  is 
presented  in  Table  V.  Statistical  confidence  in  w'r'  exists 


TABLE  V 

CU  absolute  humidity  flux  composite  and 
corresponding  standard  deviations. 


w'r'  (gs”^m~^) 
(std.  dev.) 

Left 

Center 

Right 

Downwind 

0.045 

(0.003) 

0.041 

(0.005) 

0.045 

(0.003) 

Middle 

0.039 

(0.005) 

0.068 

0.041 

(0.008) 

Upwind 

0.049 

(0.003) 

0.053 

(0.010) 

0.049 

(0.003) 
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for  all  bins  of  the  composite.  As  is  expected  by  parcel 
displacement  theory  for  a  CMASL,  flux  in  all  bins  of  the 
composite  is  down-gradient.  This  property  is  also  found  in 
the  majority  of  updrafts  and  downdrafts  of  Khalsa  and  Greenhut 
(1985) .  Composite  patterns  of  w'r'  resemble,  for  the  most 
part,  those  of  w'w'.  The  largest  flux  magnitude  is  found  in 
the  middle  center  bin.  Fluxes  along  the  upwind  side  of  the 
draft  is  slightly  larger  than  those  on  the  downwind  side. 
There  is,  however,  symmetry  along  the  cross-mean-wind  group 
axis. 

The  CD  w'r'  composite  is  presented  in  Table  VI.  As  with 


TABLE  VI 

CD  absolute  humidity  flux  composite  and 
corresponding  standard  deviations. 


w'r'  (gs~^m“^) 
(std.  dev.) 

Left 

Center 

Right 

Downwind 

0.025 

(0.003) 

0.025 

(0.003) 

0.025 

(0.003) 

Middle 

0.023 

(0.004) 

0.034 

0.022 

(0.003) 

Upwind 

0.020 

(0.002) 

0.021 

(0.002) 

0.020 

(0.002) 
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the  CU,  there  is  statistical  confidence  in  the  entire 
composite  and  the  absolute  humidity  flux  is  down-gradient 
throughout  the  composite,  with  a  maximum  at  the  center.  There 
are,  however,  some  differences  from  the  CU  w'r'  composite.  As 
is  expected  from  parcel  displacement  theory  and  the  relative 
magnitudes  of  flux  in  the  w'w'  draft  composites,  a  weaker 
absolute  humidity  flux  composite  is  observed  for  the  CD.  Also, 
the  downwind  side  of  the  CD  composite  contains  slightly 
stronger  flux  than  that  of  the  upwind  side  but  as  in  the  CU 
composite,  symmetry  exists  along  the  cross-mean-wind  group 
axis. 

3.4  MOMENTUM  FLUX 

There  is  statistical  confidence  for  all  bins  of  the  CU 
planview  along-mean-wind  momentum  flux  composite,  presented  in 
Table  VII.  The  flux  contained  in  all  bins  are  down-gradient, 
given  the  usual  sign  of  the  surface  layer  wind  shear.  This 
same  property  is  seen  in  moist  updrafts  by  Khalsa  and  Greenhut 
(1985) .  Wilczak  (1984) ,  however,  shows  a  region  of  counter¬ 
gradient  momentum  flux  just  inside  the  trailing  edge  of  his 
ensemble  updraft  vertical  cross  section.  The  largest  flux 
magnitude  in  the  CU  composite  is  found  in  the  middle  center 
bin.  Less  significant  asymmetries  exist.  For  example,  w'u' 
is  stronger  along  the  downwind  side  of  the  composite  than 
along  the  upwind  side.  Weakest  flux  is  found  in  the  middle 
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left  and  right  bins. 

The  CD  planview  flux  composite  of  w'u'  is  presented  in 
Table  VIII.  Statistical  confidence  exists  in  all  data  with 


TABLE  VII 

CU  along-mean-wind  momentum  flux  composite  and 
_ corresponding  standard  deviations. _ 


w'u'  (m^s"^) 
(std.  dev.) 

Left 

Center 

Right 

Downwind 

-0.133 

(0.008) 

-0.145 

(0.041) 

-0.133 

(0.008) 

Middle 

-0.107 

(0.016) 

-0.219 

-0.099 

(0.018) 

Upwind 

-0.122 

(0.023) 

-0.123 

(0.037) 

-0.122 

(0.023) 

the  exception  of  the  upwind  center  bin,  where  the  standard 
deviation  approaches  the  magnitude  of  the  flux.  As  with  the 
CU,  and  as  seen  in  Khalsa  and  Greenhut  (1985) ,  w'u'  is  down- 
gradient  throughout  the  composite.  As  discussed  below,  this 
finding  differs  from  that  of  Wilczak  (1984)  for  his  ensemble 
downdraft  vertical  cross  section  of  w'u'.  The  region  of 


strongest  flux  within  the  CD  composite  is  found  on  the 
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downwind  side,  while  weakest  flux  is  on  the  upwind  side.  The 
largest  magnitude  of  w'u'  is  seen  in  the  downwind  center  bin 
while  the  smallest  is  in  the  upwind  center  bin.  Cross-mean- 
wind  group  axis  symmetry  is  seen. 

As  stated  above,  Wilczak  (1984),  using  finer 
observational  resolution  than  was  possible  in  this  study, 
shows  a  region  of  counter-gradient  momentum  flux  just  inside 


TABLE  VIII 

CD  a long-mean-wind  momentum  flux  composite  and 
_  corresponding  standard  deviations. _ 


w'u'  (m^s“^) 

( std .  dev . ) 

Left 

Center 

Right 

Downwind 

-0.077 

(0.014) 

-0.092 

(0.027) 

-0.077 

(0.014) 

Middle 

-0.065 

(0.015) 

-0.084 

-0.063 

(0.016) 

Upwind 

-0.043 

(0.008) 

-0.028 

(0.022) 

-0.043 

(0.008) 

the  trailing  edge  of  his  ensemble  updrafts  and  downdrafts. 
Given  the  differences  in  horizontal  resolution,  it  is  possible 
that  the  along-mean-wind  asymmetries  in  both  studies  are 


related 
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4 .  Summary 

The  coupling  of  the  atmosphere  to  the  ocean  surface  is 
accomplished  predominately  through  fluxes  induced  by  drafts 
within  the  marine  atmospheric  surface  layer.  Planview 
composite  draft  eddy  correlation  flux  patterns  of  vertical 
velocity,  buoyancy,  absolute  humidity,  and  the  along-mean-wind 
component  of  momentum,  within  the  convective  marine 
atmospheric  surface  layer,  are  derived  in  this  study.  High 
resolution  turbulence  data  are  obtained  from  20  NCAR  Electra 
50  m  mean  sea  level  flight  legs  from  Project  FIRE.  Data 
processing  results  in  a  demeaned,  detrended,  despiked,  high- 
pass  filtered  perturbation  dataset. 

A  conditional  sampling  technique  is  used  to  isolate 
updrafts  and  downdrafts  of  interest  from  the  environment.  The 
technique  uses  horizontal  scale  and  magnitude  thresholds  based 
on  w'  event  criteria.  The  w'  dataset  is  band-pass  filtered  to 
eliminate  the  inertial  subrange  contribution  (0.1  Zj^  cutoff 
wavelength)  and  the  mesoscale  contribution  (1.5  cutoff 
wavelength)  .  The  minimum  draft  width  is  50  m  while  0.1  m/s  is 
the  w'  magnitude  threshold.  This  technique  results  in  the 
isolation  of  2839  updrafts,  '  -aging  98  m  in  width,  and  3062 
downdrafts,  averaging  107  m  in  width. 

After  the  identification  of  drafts  of  interest,  eddy 
correlation  fluxes  are  calculated  for  those  portions  of  the 
processed  dataset  where  drafts  of  interest  are  located.  Mean 
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wind  positioning  within  drafts,  draft  symmetry,  and  mean  draft 
intensity  are  used  to  implement  a  series  of  averaging 
calculations,  grouping  techniques,  and  scaling  operations. 

The  result  of  the  above  data  manipulation  is  a  planview 
composite  updraft  and  downdraft  for  each  flux  discussed.  The 
composites  are  composed  of  flux  data  located  in  horizontal 
three  bin  by  three  bin  arrays.  Each  bin  represents  a  mean- 
wind-relative  sector  of  a  composite  draft:  upwind,  middle,  or 
downwind  by  left,  center,  or  right.  Along-mean-wind  group 
axes  extend  through  the  center  of  composite  drafts  while 
cross-mean-wind  group  axes  extend  through  the  middle  of 
composite  drafts. 

Statistical  confidence  exists  in  the  two  vertical 
velocity  flux  composites.  The  strongest  flux  region  of  the 
updraft  composite  is  located  along  the  along-mean-wind  group 
axis,  with  the  largest  magnitude  flux  found  at  the  center  of 
the  draft.  While  the  largest  magnitude  of  w'w'  within  the 
downdraft  composite  is  also  found  at  the  center  of  the  draft, 
the  overall  composite  is  weaker  and  more  uniform  than  that  of 
the  updraft.  In  the  buoyancy  flux  composites,  only  data  in 
bins  along  the  along-mean-wind  and  cross-mean-wind  group  axes 
are  statistically  significant.  Along  these  axes,  symmetries 
exist  and  the  fluxes  are  down-gradient.  Statistical 
confidence  exists  throughout  both  absolute  humidity  flux 
composites.  As  is  expected  by  parcel  displacement  theory,  all 
fluxes  are  down-gradient.  The  largest  magnitudes  of  flux  are 
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found  in  the  middle  center  bins,  and  both  composites  show 
cross-mean-wind  group  axis  symmetry.  The  updraft  and 
downdraft  along-mean-wind  component  of  momentum  flux 
composites  are  statistically  significant  throughout,  with  the 
exception  of  the  upwind  center  bin  of  the  downdraft  composite. 
All  fluxes  of  both  composites  are  down-gradient.  The  largest 
flux  magnitude  for  the  updraft  is  in  the  middle  center  bin 
while  that  of  the  downdraft  is  located  in  the  downwind  center 
bin.  Along-mean-wind  asymmetry  exists  within  both  composites 
with  downwind  bins  having  greater  magnitudes  of  flux  than 
upwind  bins. 

It  is  hoped  that  this  study  provides  new  useful  insight 
into  the  planview  flux  structure  of  CMASL  updrafts  and 
downdrafts.  The  degree  of  spatial  resolution  in  these 
results,  as  well  as  the  techniques  presented  in  this  study, 
could  be  of  use  in  verification  of  nonlinear  convective 
boundary  layer  models  such  as  that  presented  in  Haack  and 
Shirer  (1992).  These  results  are  also  appropriate  for  the 
input  of  forcing  in  time  dependent  two-dimensional  ocean  wave 
models  of  cats  paw-type  features. 

Of  course  much  room  exists  for  improvement  upon  and 
expansion  beyond  this  study.  This  study  should  be  used  as  a 
stepping  stone  for  such  research  initiatives. 
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or  THI  MARINI  ATMOSPKCRIC  SURFACt  LAYER 


Todd  0.  Sikora  and  daorga  S.  Young 

Tha  Pannaylvania  StaM  onlvaralty 
Univaraity  Park,  Pannaylvanla 


1 .  INTRODOCriOM 

Alr/aaa  fluxaa  of  haat,  nolatura,  and 
fflORiantum  ara  raalizad  pradonlnantly  through  aub- 
maaoacala  convactivaly  drlvan  updrafta  (CDUs)  and 
downdrafts  (CDDs),  dianatara  of  which  ara  on  tha 
ordar  of  tans  to  thousands  of  macars  (Lansehow 
and  Stsphans,  1980;  Khalsa  and  Graanhut,  1985). 
Various  analysis  niathods  may  ba  amployad  to  study 
thasa  drafts.  Ona  approach  is  to  obsarva  tha 
affact  of  tha  drafts  on  othar  aspacta  of  tha. 
atmosphara  or  ocaan,  such  as  cat'a-paw  wava 
packats  or  cumulus  clouds.  A  mora  quantitativa 
approach,  takan  in  this  study,  involvaa  sampling 
tha  turbulant  fluctuations  of  tamparatura, 
moistura,  and  valocity  within  tha  drafts. 
Conditional  sampling  of  obsarvad  or  laodalad 
turbulanca  fialds  has  baan  proven  to  ba  affactiva 
in  tha  quantitativa  analysis  of  convactiva  draft 
structura  and  dynamics  (Lansehow  and  Stephans, 
1980;  Hilesak  and  Tillman,  1980;  Graanhut  and 
Khalsa,  1982;  Wilczak,  1984;  Khalsa  and  Graanhut, 
1985;  Young,  1988;  Schumann  and  Moang,  1991). 

Whila  this  pravioua  work  haa  raaultad  in 
datailad  vertical  profiles  of  propartias  of 
marina  atmoapharic  boundary  layac  (KABL)  updrafts 
and  downdrafts,  information  on  planviaw  flux 
variability  of  thasa  KABL  CDUs  and  cODs  is 
lacking.  Our  study  uses  conditionally  sampled 
and  eompositad  eddy  correlated  FIRE  (First  ISCCP 
(International  Satallita  Cloud  Climatology 
Project)  Regional  Sxpacimant)  aircraft  data  to 
investigate  tha  typical  cou  and  COD  planviaw  flux 
patterns  of  vertical  valocity  (w'),  buoyancy 
(t  '),  absolute  humidity  (r'),  and  along-maan-wind 
momantum  (u')  in  tha  convective  marina 
atmospheric  surface  layer  (CMASL). 

2.  PROCEDURES 

The  usa  of  NCAR  Elactra  turbulanca 
maasuramants  from  Project  FIRE  is  advantageous 
for  our  analysis.  During  FIRE,  tha  Elactra 
collactad  high  rasolutlon  turbulanca  data  from  a 
large  number  of  drafts  within  the  CMASL.  While 
each  draft  was  penetrated  only  once,  panatrations 
of  many  different  drafts  ware  made  from  numerous 
angles  relativa  to  tha  mean  wind.  Data  ware 
obtained  from  20  horizontal  lags  flown  at  50  m 
MSL.  Bach  lag  was  approximately  SO  km  (500  sac) 
in  length.  Thasa  lags  ware  flown  over  a  period 
of  seven  days  in  July,  1987  off  tha  coast  of 
California.  The  boundary  layer  depth  (Z,) 
averaged  920  m  during  these  legs  so  the  50  m 
flight  level  was  well  within  the  surface  layer. 
Ratios  of  the  boundary  layer  depths  to  tha  Monln- 
Obu)chov  lengths  ranged  from  -0.62  to  -37.56, 
indicative  of  slightly  to  moderately  unstable 
KABLs. 

This  study  uses  the  same  instrusMntation 
.reviewed  by  Nucciarone  and  Young  (1991).  Data 
preprocessing  foe  mean  and  trend  resnval 
parallels  that  of  Lansehow  and  Stephens  (1980). 
Masoscala  phenomena  are  eliminated  from  the  data 
series  using  tha  highpass  filtering  techniques  of 
Young  (1987,  1988)  incorporating  a  1.5  S,  cutoff 
wavelength  based  on  the  spectra  of  Nucciarone  and 
Young  (1991).  Coordinate  transformations  are 
perfonaad  on  the  resulting  perturbation  wind 
components  yielding  ona  in  the  direction  of  the 
mean  wind,  u',  and  another  aligned  with  positive 
values  90*  to  tha  left  of  tha  mean  wind  (v*).  The 
latter  is  both  statistically  and  physically 
insignificant  (is.  small  in  magnitude)  whan 
compared  with  tha  alongwind  component  and, 
therefore,  its  flux  is  not  discussed. 


Conditional  sampling  to  distinguish 
updrafts  and  downdrafts  uses  event  criteria  based 
on  w',  following  Graanhut  and  Khalsa  (1982),  and 
incorporates  garmane  horizontal  scala  and 
magnitude  thresholds,  following  Young  (1988). 
Thus,  for  event  identification,  tha  w'  series  is 
bandpass  filtered  to  eliminate  tha  inertial 
subrange  contribution  (0.1  Z^  cutoff  wavelength) 
and  tha  masoscala  contribution  (1.5  Z,  cutoff 
wavelength  as  mentioned  above).  Tha  inertial 
subrange  contribution  is,  in  contrast,  retained 
in  tha  data  series  used  for  eddy  correlation 
statistic  computation.  Given  tha  average  Z,  of 
this  study,  a  minimum  draft  width  of  50  m  is 
chosen  so  that  tha  ahortast  acceptable 
u^raft /downdraft  couplat  corresponds  with  tha 
0.1  Z,  inertial  subrange  cutoff  wavelength  used  in 
tha  bandpass  filter.  Inspection  of  tha  FIRE  w' 
data  series  suggests  that  a  w'  threshold  of  ♦  0.1 
m/s  aliminatas  disorganized  areas  of  weak  ascent 
or  descant.  0.1  m/s  is  tharafora  used  as  tha 
magnitude  threshold  in  our  study.  Any  data 
series  segment  not  masting  tha  above  criteria  is 
not  considarad  to  ba  part  of  a  significant 
convective  draft  and  is  rejected.  This  method 
yields  2839  CDUs  averaging  98  o  in  width  and  3062 
CDDs  averaging  107  m  in  width. 

After  tha  isolation  of  CDUs  and  CDDs  from 
each  lag,  as  discussed  abova,  eddy  correlation 
fluxes  ara  calculated  for  each  draft.  In  ordar 
to  develop  typical  planviaw  flux  snapshots,  tha 
fluxes  from  all  the  drafts  ara  composited 
together  using  a  secies  of  averaging  calculations 
and  grouping  techniques.  Bach  draft  is  first 
divided  into  three  bins  of  equal  length;  the 
aircraft  entry  region  of  the  draft,  tha  middle  of 
the  draft,  and  the' aircraft  exit  region  of  the 
draft.  The  average  CDU  and  CDD  flux  patterns  for 
each  leg  are  then  found  by  calculating  the  bln 
average  eddy  correlation  statistics. 

The  resulting  40  leg  average  drafts  from 
the  20  flight  legs  are  then  separated  into  three 
equal  angular-width  groupe  based  on  the  path 
taken  by  the  aircraft  through  the  draft,  relative 
to  the  mean  wind.  In  polar  coordinated  aligned 
with  the  mean  wind  blowing  towards  0  ,  these 
relative  heading  groups  are;  alongwind  (330  to 
30*),  crosswind  (60*  to  120“),  and  diagonal  to  the 
wind  (30*  to  60*).  To  maintain  relative-wind 
consistency,  the  entry  and  exit  bins  are  swapped 
for  headings  opposite  of  those  given  above. 
Draft  symmetry  around  the  alongwind  axis 
(Wilczak,  1984)  justifies  combining  the  data  from 
all  four  sectors  diagonal  to  the  wind  into  the 
one  given  above.  Our  crosswind  composites  verify 
this  symmetry.  Corresponding  bln  averages  for 
each  heading  group  are  composited,  yielding  a 
typical  CDU  flux  pattern  and  a  corresponding  COO 
flux  pattern  for  each  flux  statistic.  Because  of 
the  need  to  correct  for  mean  draft  intensity 
differences  between  flight  lege,  the  perimeter 
bin  fluxes  are  rescaled  before  compositing.  This 
resesling  uses  a  correction  coefficient  based  on 
the  requirement  for  similarity  in  the  middle 
bins,  where  all  heading  groups  are  assumed  to 
overlap. 

3.  RESULTS 

Planviaw  draft  flux  compoj(^M  of  vertical 
velocity  (w*w*),  buoyancy  (w‘t^')»  absolute 
humidity  ('w*?') ,  and  along-mean-wind  momentum 
(JTiT),  are  found  in  Tables  1  through  4.  The 
tables  labeled  (A)  contain  COO  flux  composites 
while  tables  labeled  (B)  contain  COD  flux 
_ _  composites.  Each  table  is  a  three  by  three 
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horizontal  array  of  fluxaa  obtainad  from  tha  bln 
avaraglng  and  grouping  tachnlquas  daacrlbad  In 
tha  pravloua  aactlon.  For  all  tablas,  tha  maan 
wind  la  dlractad  towards  tha  top  of  tha  paga. 

Tha  planvlaw  CDU  and  CDD  vartlcal  valoclty 
flux  conpoaltas  ara  shovm  In  Tablas  lA  and  1B> 
Bacausa  all  bln  standard  davlatlona  ara  much 
amallar  than  tha  corraspondlng  fluxas  (lass  than 
1/3  of  thasa  aiaans ) .  thara  la  statistical 
confldanca  In  tha  vTvT  compoaltas.  Both  flux 
compoaltas  ara  naarly  synmatclc  across  tha  axis 
of  tha  maan  wind.  Tha  strongast  flux  raglon  In 
tha  CDU  composlta  Is  found  on  tha  alongwlnd  axis 
with  tha  largast  magnltuda  of  flux  balng  at  tha 
draft  cantar.  This  shapa  comparas  wall  with  that 
found  at  SO  m  ovar  land  by  Wllczak  (1984).  whlla 
tha  largast  magnltuda  of  w‘w‘  In  tha  CDO 
composlta  Is  also  found  In  tha  cantar  bln,  tha 
ovarall  flux  composlta  la  mora  uniform  than  that 
of  tha  CDU  and  la  waakar.  In  agraamant  with 
findings  of  Graanhut  and  Khalsa  (19B2). 

Tabla  lA.  CDU  vartlcal  valoclty  flux  composlta. 


w'w'  (m*s'*) 

Left 

Cantar 

Right 

Downwind 

0.227 

0.244 

0.227 

Middle 

0.217 

0.341 

0.225 

Upwind 

0.238 

0.266 

0.238 

Tabla  IB.  COD  vartlcal  valoclty  flux  composlta. 


1  w'w'  (m*s‘*) 

Laft 

Cantar 

Right 

1  Downwind 

0.160 

0.161 

0.160 

1  Mlddla 

0.171 

0.204 

0,173 

1  Upwind 

0.162 

0.166 

0.162 

Tha  buoyancy  flux  of  tha  currant  study 
approxlmatss,  to  a  larga  axtant .  tha  hast  flux. 
Tha  planvlaw  CDU  and  COD  w'ty'  compoaltas  arm 
shown  In  Tablas  2A  and  2B.  Thara  la  statistical 
confldanca  In  thosa  parts  of  tha  compoaltas  mada 
up  of  tha  along  wind  and  crosswind  groups  but  not 
In  tha  diagonal  alamants  (array  cornars). 
Alongwlnd  and  crosswind  symmatrlas  axlst  for  both 
tha  CDU  and  COD  composlta.  All  statistically 
significant  fluxas  within  tha  compoaltas  ara 
down-gradlant.  In  agraamant  with  tha  lowaat 
lavals  of  Khalsa  and  Graanhut  (1985). 


(1985).  Similar  to  tha  ww'  compoaltas,  tha 
largast  flux  magnltuda  Is  found  In  tha  mlddla 
bln.  Flux  along  tha  upwind  slda  of  the  CDU 
composlta  and  downwind  slda  of  CDD  composlta  Is 
slightly  largar  than  that  of  tha  oppoalta  slda  of 
aach  raspactlva  composlta.  Thara  Is,  howavar, 
good  crosswind  symmatry  In  both. 

Tabla  3A.  CDU  absoluta  humidity  flux  composite. 


”'^'l  -2 
(gs  ’m  *) 

Left 

Cantar 

Right  1 

Dovmwlnd 

0.045 

0.041 

0.045 

1  Middle 

0.039 

0.068 

0.041 

8  upwind 

0.049 

0.053 

0.049  ’ 

Tabla  3B.  CDD  absoluta  humidity  flux  composlta. 


Cantar 

Right  1 

1  Downwind 

0.025 

0.025 

0.025 

1  Mlddla 

0.023 

0.034 

0.022 

1  Upwind 

0.020 

0.021 

0.020 

Tha  CDU  and  CDD  planvlaw  w'u'  compoaltas 
ara  shown  In  Tablas  4A  and  4B.  Thara  Is 
statistical  confldanca  In  tha  fluxas  of  all  bins 
axcapt  that  of  tha  upwind  cantar  bln  In  tha  CDD 
composlta.  All  fluxas  of  both  composites  ara 
down-gradlant,  given  tha  usual  surface  layer 
vartlcal  wind  shear.  This  same  property  Is  seen 
In  similar  drafts  of  Khalsa  and  Graanhut  (1985). 
Tha  largest  CDD  flux  magnltuda  Is  again  at  tha 
canter  of  tha  draft  whlla  tha  CDD  flux  magnltuda 
maximum  is  located  at  tha  downwind  center.  Tha 
CDU  flux  magnitude  does,  howavar,  show  soma 
tendency  towards  this  form  of  asymmetry  with  tha 
downwind  bins  having  graatar  magnitudes  than  tha 
upwind  bins.  Mllczak  (1984),  using  finer 
horizontal  resolution,  shows  a  raglon  of  counter- 
gradient  momentum  flux  just  Inside  tha  upwind 
edge  of  his  ansambla  updrsfts  and  downdrafts. 
Civan  tha  differences  In  horizontal  resolution. 
It  Is  possible  that  tha  alongwlnd  asymmetries  In 
both  studies  ara  related. 

Tabla  4A.  CDU  along-maan-wlnd  momentum  flux 
com]poslta. 


Tabla  2A.  CDD  buoyancy  flux  composlta. 


Laft 

Cantar 

Right 

Downwind 

0.024 

0.010 

0.024 

Mlddla 

0.010 

0.018 

0.011 

Upwind 

0.006 

0.014 

0.006 

Tabla  28.  CDO  buoyancy  flux  composlta. 


mm 

Laft 

Cantar 

Right 

Downwind 

0.009 

0.005 

0.009 

Middle 

0.004 

0.006 

0.004 

Upwind 

-0.002 

0.004 

-0.002 

Laft 

Cantar 

Right 

Downwind 

-0.133 

-0.145 

-0.133 

Mlddla 

-0.107 

-0.219 

-0.099 

Upwind 

-0.122 

-0.122 

Tabla  4B.  COO  along-BWan-wlnd  momantum  flux 
composlta. 


■cieJiflUB 

Laft 

Cantar 

■m 

Downwind 

-0.077 

-0.092 

-0.077 

Middle 

-0.084 

-0.063 

Upwind 

-0.043 

-0.028 

-0.043 

Tha  planvlaw  CDU  and  CDO  absoluta  humidity 
flux  cos^sitas  ara  shown  In  Tablas  3A  and  38. 
Statistic^  confldanca  of  w^r'_axlsts  in  all  bins 
of  both  compoaltas.  As'  Is  sxpact^'  by  parcal 
dlsplacamant  thaory,  tha  flux  composites  ara 
down-gradlant,  as  Is  found  In  tha  sMjorlty  of 
updrafts , and  downdrafts  of  Khalsa  and  Graanhut 


4 .  ACKMOMLBDGBMKirrS 

Ha  wish  to  thank  Ors.  Khalsa  and  Shlrar  for 
valuable  discussions.  Funding  for  this  research 
has  baan  through  tha  Office  of  Haval  Research 
undar  grants  N00014-90-J-4012  and  N00014-91-J- 
1992. 


26 


S .  MFERINCXS 

Graanhut,  C.  K. ,  «nd  S.  J.  S.  Khalaa,  1982i 
Updraft  and  downdraft  avanta  in  tha 
atinoapharic  boundary  layar  ovar  tha 
aquatorlal  Pacific  Ocaan.  J.  Aemoa.  Sci.,  39, 
1803-1818. 

Khalsa,  S.  J.  S.,  and  C.  K.  Graanhut,  1985 t 
Conditional  aanpling  of  updrafta  and 
downdrafts  in  tha  marina  atmoapharlc  boundary 
layar.  J.  Atmos.  Sci.,  43,  2550-2562. 

Lanschow,  D.  H.,  and  P.  L.  Staphana,  1980t  Tha 
rola  of  tharmals  in  tha  convactiva  boundary 
layar.  Bound-L4y0r  Macaor.,  19,  509-532. 

Nucciarona,  J.  J. ,  and  G.  S.  foung,  1991: 

Aircraft  maasuramanta  of  turbulanca  spactra 
in  tha  marina  atratocumulua-toppad  boundary 
layar.  J.  Atmos.  Sci.,  48,  2383-2393. 

Schumann,  U.  and  C.-H.  Moang,  1991:  Pluma  fluxas 
in  claar  and  cloudy  convactiva  boundary 
layara.  J.  Atmos.  Sci.,  48,  1746-1757. 

Wilczak,  J.  M.,  1984:  Larga-scala  addias  in  tha 
unstably  stratifiad  atmoapharlc  surfaca 
layar.  Part  I:  Valocity  and  tamparatura 
atiructura.  J.  Atmos.  Sci.,  41,  3537-3550. 

Wilczak,  J.  M.,  and  J.  X.  Tillman,  1980:  Tha 

thraa-dimanaional  atructura  of  convaction  in 
tha  atmoapharlc  aurfaca  layar.  J.  Atmos. 
Sci.,  37,  2424-2443. 

young,  G.  S.,  1987:  Hixad  layar  spactra  from 
aircraft  maasuramanta.  J.  Atmos.  Sci.,  44, 
1251-1256. 

Young,  a.  S.,  1988:  Turbulanca  atructura  of  tha 
convactiva  boundary  layar.  Part  II:  Phoanix 
78  aircraft  obsarvationa  of  tharmals  and 
thair  anvironmant.  J.  Atmos.  Sci.,  45,  727- 
735. 


27 


OFFICE  OF  NAVAL  RESEARCH 

PUBLI CATIONS / PATENTS / PRES ENTATIONS / HONORS  REPORT 
1  Oct;  91  through  3  0  Sep  92 


R&T  Number; 


Contract/Grant  Number: 
Contract/ Grant  Title: 
Principal  Investigator: 
Mailing  Address: 


4215021 - 03 

N00014-90-J-4012 

SAR-Related  Stress  Variability  in  the  Marine 
Atmospheric  Boundary  Layer  (MABL) 

H.N.  Shirer  and  G.S.  Young 

503  Walker  Building,  Department  of  Meteorology, 
Penn  State  University,  University  Park,  PA  16802 


Phone  Number:  (with  Area  Code)  814-863-1992  (HNS);  814-863-4228  (GSY) 


E-Mail  Address: 


HNSiapSUVM. PSU.EDU;  G.Young/Omnet 


a.  Number  of  Papers  Submitted  to  Refereed  Journal  but  not  yet 

published:  0 

b.  Number  of  Papers  Published  in  Refereed  Journals:  ^ 

(list  attached) 

c.  Number  of  Boolcs  or  Chaoters  Submitted  but  not  vet  Ptiblished: 


d.  Number  of  Books  or  Chaoters  Published; 


(list  attached) 


e.  Number  of  Printed  Technical  Reports  &  Non-Ref ereed  Papers: 
(list  attached) 


Number  of  Patents  Filed: 
Nximber  of  Patents  Granted: 


(list  attached) 


h.  Number  of  Invited  Presentations  at  Worlcshops  or  Prof.  Society 

Meetings :  6 

.  2 

i.  Number  of  Presentations  at  Workshops  or  Prof.  Society  Meetings:  _ 

j.  Honors /Awards /Prizes  for  Contract/ Grant  Employees:  none 
(list  attached,  this  might  include  Scientific  Soc.  Awards/ 

Offices,  Pr  emotions  /  Faculty  .Aw2u:ds/Off  ices,  etc.) 

k.  Total  number  of  Graduate  Students  and  Post-Docs  Supported  at 
least  25%  this  year  on  this  contract /grant: 

Grad  Students  ^  and  Post-Docs  ^  including 
Grad  Student  Female  Q  and  Post-Docs  Female  ^  _ 

Grad  Student  Minority  1  and  Post-Doc  Minority  0 _ 

Minorities  include  Blacks,  Aleuts,  Amindians,  Hispanics,  etc. 

N3:  Asians  are  not  considered  an  under-represented  or  minority  group  in 

science  and  engineering. 


Enclosure  (3) 


Refereed  Papers 


Sikora,  T.D.  and  G.S.  Young,  1993:  Observations  of  planview  flux  patterns 
within  convective  structures  of  the  marine  atmospheric  surface  layer. 
Boundary  Layer  Meteorology,  35,  in  press. 


Non-Refereed  Papers 

Hare,  J.E,  1992:  Shipboard  eddy-covariance  measurements  of  the  turbulent 
fluxes  of  heat,  moisture,  and  momentum.  MS  Thesis,  Penn  State 
University,  207  pp. 

Sikora,  T.D.,  1992:  Air/sea  flux  patterns  within  convective  structures  of  the 
marine  atmospheric  surface  layer,  MS  Thesis,  Penn  State  University,  28 

pp. 

Sikora,  T.D  and  G.S.  Young,  1992,  Air/sea  flux  patterns  within  convective 
structures  of  the  marine  atmospheric  surface  layer.  Preprints,  Tenth 
Symposium  on  Turbulence  and  Division,  Sept.  29-Oct.  2,  1992,  Portland, 
OR,  American  Meteorological  Society,  25-27. 


Planview  ©f  Composite 
Downdraft  Horizontal  Wind 


Planview  of  Composite 
Updraft  Horizontal  Wind 


Per^urbatiomi  Wind  Fi' 


